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Abstract 

Qn 1 We study the structure of vortex solutions in a Ginzburg-Landau system for two 

■ complex valued order parameters. We consider the Dirichlet problem in the disk in R 2 
with symmetric, degree-one boundary condition, as well as the associated degree-one 
entire solutions in all of R 2 . Each problem has degree-one equivariant solutions with 
radially symmetric profile vanishing at the origin, of the same form as the unique (com- 
plex scalar) Ginzburg-Landau minimizer. We find that there is a range of parameters 
for which these equivariant solutions are the unique locally energy minimizing solu- 

■ tions for the coupled system. Surprisingly, there is also a parameter regime in which 
the equivariant solutions are unstable, and minimizers must vanish separately in each 
component of the order parameter. 
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1 Introduction 

We continue our study of the structure of vortices in two-component Ginzburg-Landau 
functionals begun in |AB2| |A"BMJ. Let Q C I 2 be a smooth, bounded domain, and ^ G 
C 2 ). We define an energy functional, 

M*; n) = I {^|v*| 2 + A (|*| 2 - 1) 2 + A (|^ + | 2 _ |^| 2 ) 2 } dx , (i) 

where ^ = (3 > and e > are parameters. Energy functionals of a form similar 

to E e have been introduced in physical models, and at the end of the section we will briefly 
describe two such contexts: a Spinor Ginzburg-Landau functional, describing ferromagnetic 
and antiferromagnetic superconductors, and giving rise to half-integer degree vortices; and 
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a Gross-Pitaevskii functional for a two- component Bose-Einstein condensate. Although the 
physical models are more complex, we expect that the essential features of the singular limit 
e — » in the physical systems will be well described by the simpler energy ([I]). 

As a model problem, we consider (pQ) with appropriate Dirichlet boundary conditions, 
and study the behavior of energy minimizers as 6 0. In the limit, minimizers \1/ should 



lie on the manifold in C 2 on which the potential F(\I>) = \(\^\ 2 - l) 2 + § (\ip+\ 2 - |^-| 2 ) 2 
vanishes. That manifold is a 2-torus E C § 3 C C 2 , parametrized by two real phases 



and thus a E- valued map ^(x) carries a pair of integer- valued degrees around any closed 
curve C, 



If the given Dirichlet boundary condition has nonzero degree in either component, then there 
is no finite energy map ^ which takes values in E and satisfies those boundary conditions, 
and we expect that vortices will be created in the e — > limit, just as in the classical 
Ginzburg-Landau model [BBHJ. 

An analysis of the global minimizers of the Dirichlet problem, with E-valued boundary 
data, is given in |AB1| IAB2] . As in the seminal work of Bethuel, Brezis, & Helein [BBHJ , 
nonzero degree boundary data give rise to vortices in Q, each of degree one in one (or 
both) of the two phases a±, and the location of the vortices is determined by minimizing a 
renormalized energy, which is derived by sharp estimates of the interaction energy between 
the vortices. The essential difference between the classical Ginzburg-Landau model and the 
energy ([T]) is that there are different species of vortices, allowing for winding in one or both of 
the two phases, a±. In the renormalized energy expansion, vortices with winding in different 
components do not interact directly, in the sense that there is no term in the renormalized 
energy which couples the location of the a + and a_ vortices. However, we discovered that 
there is a very short-range interaction between these two species due to the energy of the 
vortex cores. In particular, it is shown (for certain values of the parameter /?,) that it may 
be beneficial for two vortices of different type (one with degree [n+, n_] = [1, 0] and one with 
[n + , n_] = [0, 1]) to coincide in the e — > limit, rather than to converge to distinct points (as 
the renormalized energy would normally dictate.) In this paper, we study the finer structure 
of these compound vortices: for e > small we ask, do they resemble Ginzburg-Landau 
vortices, with = at a common vortex location; or does each component ip± vanish 
separately? 




deg(^;C) = [n+,7i_] 



n+ = deg(^+; C) 



n_ = deg(^_; C). 
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To illustrate, we begin with a simple but prototypical example: let = Di = D(0, 1), 
the unit disk, and denote 

We consider minima (or more generally, critical points) of E e over the space H, consisting 
of all functions $ G C 2 ) with the symmetric boundary condition: 

= ^ (e'V*) • (2) 

The degree on the boundary is [n + ,n_] = [1, 1], and the energy expansion of |AB2] shows 
that the minimizers \l/ e converge to a S-valued harmonic map in D x with a single limiting 
vortex at the origin. Thus, minimizers produce a compound vortex with winding in each 
phase near the origin. For e > small but nonzero, do the zeros of the two components 
coincide or not? It is easy to verify that if u e minimizes the classical Ginzburg-Landau 
energy with symmetric boundary condition, 

G e (u) = [ 

then U e = 77g(^e, w e ) is a critical point of E e with boundary data (j2J). Is the Ginzburg-Landau 
minimizer U e minimizing for E t l 
We prove the following: 

Theorem 1.1. Let u e minimize G e with u £ \qo 1 = e id and 

(i) If j3 > 1, then U e minimizes E e with Dirichlet condition (j2J) for every e > 0. 

(ii) If < j3 < 1, then for all sufficiently small e > 0, U e is not the minimizer of E t with 
boundary condition (J2J). 

The proof of Theorem 11.11 is based on comparisons between the vortex core energies (see 
( TT3]) ). and on our previous results in |AB1} IABM] . and is the content of Section [2j For 
< (3 < 1, the form of minimizers is more complex and more interesting. Our results (see 
Corollary 12.51 and Lemma |2.6[ ) suggest that the compound [n+, n_] = [1,1] vortex which 
appears in the e — > limit actually breaks down into two distinct simple vortices for e > 0. 

Although this example seems quite special, in fact the symmetric minimization problem 
plays an important role in the expansion of the energy in the method of jBBHj . and we 



I|v«|" + -L(M»-i)> 
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expect that minimizers of this problem accurately describe the structure of vortices near the 
vortex core. 

Another approach to the core structure of vortices is obtained by blowing up the solution 
at scale e around the vortex center. After rescaling and passing to the limit, one obtains an 
entire solution in all of R 2 to the elliptic system, 

-A^ + = (1 - \^ + + /3(|^_| 2 - |^ + | 2 )^+,1 
-A^_ = (1 - |*| 2 )V>- - /3(|^| 2 - |^+| 2 )^.J 

Solutions to ([3]) obtained by blowing up minimizers in Q will satisfy an integrability condition, 

f \(\^-l) 2 + ^(\^-\^) 2 \dx<oo, (4) 

analogous to the condition of [BMR] for the classical Ginzburg-Landau equation. While they 
have infinite energy measured in the whole of M 2 , they do inherit a local energy minimizing 
property, identified by De Giorgi. For \1/ e if/ oc (M 2 ; C 2 ) satisfying (@J, we denote 

£(*; n) := E^; O) = J ji|W| 2 + \ (|*| 2 - l) 2 + | (|^_| 2 - |^ + | 2 ) 2 } . (5) 

Definition 1.2. We say that ^ is a locally minimizing solution of (Tj|) if holds and if for 
every bounded regular domain Q C M 2 , 

holds for every $ = (<p + , (p_) G H 1 ^; C 2 ) with &\ dn = ^f\ 9n . 
For the classical Ginzburg-Landau equation in M 2 , 

—A-u = (1 — |m| 2 )m, 

the locally minimizing solutions are completely known. Combining results by Shafrir [Sh], 
Sandier [Sa], and Mironescu |M2j . the unique nontrivial locally minimizing solution is (up 
to symmetries) the degree-one equivariant solution, u = f{r)e %e . 

In |ABM] we proved several results on the entire solutions of ([3]). Following the work of 
[BMR] on the Ginzburg-Landau equations, any solution of (jBD satisfying (BJ has a degree 
pair at infinity, (see |ABM] ) . n± = deg('0±; oo) = deg (jjj^] Sr \ for all sufficiently large radii 
R. We also showed that there exists a unique equivariant solution to (j3J) for each degree 
pair [n+, n_], but as in the Ginzburg-Landau case we do not expect all those solutions to be 
local minimizers. Indeed, it is only for the simplest, "ground state" degrees [n+,n_] = [1,0] 
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or [0, 1] that we can assert the existence of locally minimizing entire solutions. In |ABM] we 
show that these vortex solutions are "coreless", that is, |^| is bounded away from zero in 
R 2 . 

For degrees [1, 1] at infinity, it is not clear whether or not a locally minimizing solution 
exists for ([3]). For instance, it is easy to verify that, if u = f(r)e ld is the symmetric, degree 
one solution of the Ginzburg-Landau equations in R 2 , then U = {-j^ u i 75^) solves (J3J) with 
degrees [n + , n_] = [1,1]. We then ask: is U a locally minimizing solution, and are there any 
others? 

We prove the following result: 

Theorem 1.3. (%) For (3 > 1, ty* is a locally minimizing solution with degree deg(^>±, 00) = 
1 if and only if 

V = -j=(u(x- a)e i<t>+ ,u(x- a)e^- ) , (6) 

where 6+ are real constants, a G R 2 is constant, and u(x) is the (unique) equivariant 
solution to the Ginzburg-Landau equation in M 2 with deg(u, 00) = 1. 

(ii) If0<(3<l and ^* is a locally minimizing solution with degree deg(ip±, 00) = 1, then 
is bounded away from zero in M 2 . Moreover, 

[ (|^ + | 2 -hA_| 2 ) 2 >7T. (7) 

In particular, for < < 1, the Ginzburg-Landau solution is not locally minimizing 
for ([3]). This implies (see Proposition 14. ip that for < /3 < 1, local minimizers with degree 
pair [n+,n_] = [1,1] must have distinct zeros in each component. In this way, a locally 
minimizing solution for G (0, 1) should resemble a gluing together of two simple vortex 
solutions (of degrees [1,0] and [0, 1]), studied in |ABM] . This has an important implication 
for the Dirichlet problem: let Q C M 2 be a bounded, smooth domain, and g : dQ — > E a 
given smooth boundary condition. We then conclude: 

Theorem 1.4. Let < /3 < 1, and suppose \l/ e minimizes E e (^f; Q) with Dirichlet boundary 
condition ^ e \dn = g- Then there is some c > such that, for all e > sufficiently small, 
|* e | >cmVL. 

In particular, the minimizer in the disk Vt — Dj with symmetric boundary condition (T5]) 
has a single Ginzburg-Landau type vortex (with both components vanishing at the origin) 
for > 1, but for < < 1 each component vanishes separately, and \^ e \ is bounded 
away from zero. From the analysis of the renormalized energy done in [AB1, AB2J, as 
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e — > the zeros of ip + and ip- must tend to the origin. It is an interesting open question 
to determine the rate at which they coalesce as e — > 0. If the mutual distance between the 
zeros in each component is of the order of e, then blowing up at scale e produces a locally 
minimizing solution to with degree pair [1, 1] at infinity. Necessarily, this local minimizer 
is non-equivariant, with separated zeros in each component. On the other hand, if no locally 
minimizing solution exists with degree pair [n + ,n_] = [1, 1], then the distance between the 
two vortices in the boundary-value problem must necessarily be much larger than e, and the 
compound vortex breaks down into a distinct pair of = [1,0] and [n + ,n_] = [0,1] 

vortices for e > 0. 

The analysis of locally minimizing solutions to fl3]) is done in sections |3] and HI It relies 
on a priori estimates of solutions in the spirit of Brezis, Merle, & Riviere [BMRJ and Shafrir 

Finally, we return to the symmetric boundary value problem ([2]). Recall that Theorem ll.il 
states that for < (3 < 1 the symmetric solutions U e = 775(^6, u e) cannot be local minimizers 
for small e. Using a bifurcation analysis, we provide a more detailed description of how the 
symmetric solutions lose stability, and the structure of the solutions near the critical value: 

Theorem 1.5. LetO < f3 <1. 

(i) There exists e@ > for which U e is a strict local minimizer of E e for e > ep, and U e is 
unstable if e < ep. 

(ii) e = €p is a point of bifurcation for critical points of E t with Dirichlet condition §2§. 
More precisely, there exists S > and a real-analytic family {(^'^ , e(t)}\ t \<s,^es 1 of 
non-equivariant solutions bifurcating from U e „ at t = 0, and these are the only non- 
equivariant solutions in a 5 -neighborhood of (U e „,ep) in H x (0,oo). Each component 
of the non-equivariant solutions $ = (ip +) ipJ) has exactly one zero if)±(z±) = 0, and 
their zeros are antipodal z_ = —z + ^ 0. 

A more detailed description of the bifurcation from the equivariant solutions U e is given 
in Theorem 15.231 Indeed, the analysis of the linearization around U e follows the same steps 
as for the degree d>2 case for the Ginzburg-Landau functional (see Mironescu |Mlj ). and 
we show that, apart from the S 1 -symmetry of the problem, the ground state eigenspace is 
simple. 

Finally, we briefly discuss two physical contexts for our results on compound vortices. 
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Fractional degree vortices 

Our original motivation for studying the functional ([T]) comes from Spinor Ginzburg-Landau 
functionals introduced in models of ferromagnetic and antiferromagnetic superconductors 
[KR ] or Bose-Einstein Condensates (BEC) [IsM] . 

Let Q C M 2 be a smooth, bounded domain, and \l/ G H l (Q;C 2 ). We define an energy 
functional, 

= ^ / n {|V*| 2 + ^ (|*| 2 - l) 2 + ^ x ^ 2 ) 2 } 
where \I/ = ■0 2 ), ^1 x ^2 = I m (^1^2), /? > and e > are parameters. The quantity 

S = ipx x ip 2 = Im {-01 ^2} 

is interpreted as the ^-component of a spin vector, which in this two-dimensional model is 
assumed to be orthogonal to the plane of £1 

As e — > 0, energy minimizers should converge pointwise to the manifold on which the 
potential term F(^f) = (\^\ 2 — l) 2 + f (/0i x ip 2 ) 2 vanishes. Since (3 > 0, we obtain a two- 
dimensional surface (a 2-torus) E C § 3 C C 2 parametrized by two real phases, <f>,u: 

E : * = G(0,w) := (e^cosw, e^sinw). 

Notice that G is doubly-periodic with minimal period G(0 + it, u ± 7r) = G(<p, u), with each 
phase executing a half cycle. For a smooth function \I/(;r) taking values in E and a simple 
closed curve C contained in the domain of \I/ we may therefore define a pair of half-integer 
valued degrees (d^, dj) corresponding to the winding numbers of the two phases around 
E. From the above observation, these degrees satisfy d^^d^ G |Z, and d<p + d^ G Z. The 
singularities which appear in energy minimizers as e — » will thus be half-integer quantized, 
giving rise to fractional degree vortices. 

The connection between the energies E e and £ e is direct: by a unitary transformation in 
the range, 

V>± := 7/= (^l ±#2), 

the two are seen to be equal, 8 e (tpi,ip2) = E e (ip + ,ip_). In these new coordinates, the frac- 
tional degree vortex for (^1,^2) with degree pair (d^,d u ) becomes an integer quantized 
vortex for with degree pair [n + ,n_] = [dj, + d u ,d<p — dj\. In particular, the min- 

imal energy fractional degree vortices with (d^dS) = (|,±|) are associated to the integer 
degrees [n + ,n_] = [1,0] and [0,1], and the Ginzburg-Landau-like vortex of degree pair 
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dw) = (1 ; 0) becomes the compound vortex [n+,n_] = [1, 1] in the new coordinates for 
^ . We note that in the new coordinates, the spin 

1 



S = ^x^ 2 = ^_\ 2 -\^ + \ 2 ) 



remains an important quantity. 



Two- component BEC 



The functional (OQ) may also be derived from the Gross-Pitaevsky energy for a rotating 
mixture of two BEC, introduced in [SC ETCj . In this model, we consider the pair \1> = 
[y>_i_,'0_] : IR 2 — > C 2 which minimizes 

= J [\\W\ 2 - Lo(*,id e *) + lv(xm 2 + i [a + |^ + | 4 + a_|^_| 4 + 26|^ + | 2 |^-| 2 ] } , 

(8) 

where the constant u represents the angular speed of rotation, V(x) the trapping potential, 
and a + , a_ > and & 6 1 material constants. Here, and throughout the paper, we use angle 
brackets to denote the real scalar product on C or C 2 , 

= Re {$ ■ ^} , $,$GC 2 ; (<p,ip) = Re {(pip}, <p,?peC. 

The energy is to be minimized over the constraints, 

\ip+\ 2 — m +i I \ip-\ 2 = m_. 



For simplicity we replace the trapping potential V by a "flat trap" in the domain Q d IR 2 
(see |CYj for example), in other words we set V(x) = but impose a (Neumann) boundary 
condition via \& G i/ 1 (fi; C 2 ). We also assume that a + = a_ =: a > and assume the two 
species are balanced, 

/|^ + | 2 = ± = /liM 2 - (9) 

Given the constraints on the L 2 -norms, we may then complete the square in the quartic terms 
of the potential by adding in constant multiples of \ip±\ 2 without changing the minimizers, 
to arrive at an energy which more closely resembles (II]), 

sm = J n {^iv*i 2 - u,{%id&) + ^m 2 -i) 2 + | o^-i 2 - i^ + i 2 ] 2 } , 

with e 2 = 2{ f +a ^ and (5 = —2 . Thus we recover the form of (pQ) with f3 > provided that 
—a < b < a. As for the single- component BEC energy (see Ignat & Millot [IMlJ ,) we expect 
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that for oo = OQ lne|), minimizers will have vortices in Q, and blow-up around vortex centers 
will result in locally minimizing entire solutions of our system fl3]). In particular, the small- 
scale structure of vortices in the two-component BEC for —a < b < a will be determined by 
our analysis of the blow-up problem ([3]). Since the parameter regime —a < b < corresponds 
to /3 > 1, the degree [n+,n_] = [1, 1] vortices will be radially symmetric, and well-described 
by the classical Ginzburg-Landau model. However, the range < b < a corresponds to 
G (0, 1), and in this case the [n + ,n_] = [1, 1] vortices will decompose into two separate 
vortex cores, with bounded away from zero in the core, the "coreless" vortices. 

2 The Symmetric Dirichlet Problems 

We begin with the fundamental boundary value problems, with symmetric data given on the 
boundary of the unit disk Q — Di: for n± G {0, 1}, define 

I[n + ,n_]M) = min G H\B i; C 2 ), #| 9Dl = -^=(e m+d , e m - e )\ , (10) 



For future use, we also define analogous values for the disk 3s centered at the origin and of 
radius 6: 

J [n+ , n -](e,5;P) := min |s e (*) : * G H^C 2 ), *\ ms = -^=(e in + e , e in - e )} (11) 

= I[n + ,n-] (j-'i P ^ 



by scaling. For comparison purposes, we also define the analogous quantity for the Ginzburg- 
Landau functional, 

I GL (e) = min {G e (u) : u G H 1 ^; C), u\ dI$1 = e ld ) , (12) 

Ge{u)= j o Q|V M | 2 + ^(| M | 2 -1) 2 ^X, 

and Jgl(c, S) = Igl(^/S) in analogy with ( fill . Our analysis depends on comparisons between 
the following vortex core energies, (see [AB2j . jBBHj). 

Q[i,i) = Jim (/[i,i](e; 0) - vr| lne|) , 

Q[i,o] = Q[o,i] = Jim ^J[i, ](e; /?) - - | lne| 
Q GL =lim {Igl{^) - tt| lne|) 

e-S>0 

From [AB2j (using the renormalized energy for the problem /ri ) ii(e; /?)) we have 

Q[i,i](P)<2Q M (J3), for all > 0. (14) 



(13) 
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We will constantly use the following elementary but useful identity: 



+ 4(/3- 1)S 2 



and we recall that the Spin is given by 



(15) 



s=~(IV>-l 2 -hM 2 )- 

Lemma 2.1. For oil /3 > 0, %,!] (e; 0) < / GL (e). For (3 > 1, J M (e;/9) = / GL (e). 
For /3 > 1 and for any minimizer \l/ e o/ J[i ) ii(e; /?), we /ictfe 

where u(x) is a minimizer for the problem Igl{^)- 

For (3 = 1 and for any minimizer \l/ e of Jfx i i](e; 1), we have 

= ^= K(^),« 2 (x)) , 

where u 1 , u 2 are minimizers for the problem Igl{^)- 

Proof. Let u e be the minimizer of Igl(c), and £/ e = -^(u e ,u e ). Then, U e is admissible for 
/[i i] (e; j5) and has spin zero, and therefore 

I {1>1] (e;P)<E e (U e ) = G e (u e ) = I GL (e), 

for each e > and > 0. 

If /3 > 1, let \l/ e = (ijj^_,ip € _) minimize I[i,i](e; and set = v20±- Using ( TT5|) . we 

have, with^ = ^(|^| 2 -|d 2 ), 

J GL (e) <mm{GM),G e (vl)} 



<I[G e «)+G e 



F e (^ e ) 



/3-1 



0-1 



S 
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Thus, i[i,i](e;/3) = I GL (e) for all e > and > 1. For /3 > 1, S £ (x) = 0, and £ e (^ e ) = 
\[G e {v e + ) + G e {v e _)]i and we conclude that each component must be a minimizer for Igl{c)- 
Since = |t>l|, the zeros of coincide. By Theorem 9.1 of Pacard and Riviere [PRJ, 
v e + = vt- 

For = 1, E e decomposes into a Ginzburg-Landau energy for each component i/j±, and 
the conclusion is immediate. □ 

For > 1, there is a net energy saving in replacing two simple vortices, of degree pairs 
[0, 1] and [1, 0], by a single compound vortex at the same limiting location: 

Lemma 2.2. For all > 1, Q[i,i](0) = Q GL < 2Q M (p). 

Proof. The assertion Qni](0) = Qgl for all > 1 follows trivially from Lemma 12.11 and 
so Q[i,o}(0) > \Qgl for > 1 follows from flU}. Note also that when 6 = 1 by (TIoT) the 
components decouple, 

J5 e (¥) = ± (G e (V2^ + ) + G e (V2^-)) , 

and hence I[i,o](e; 1) = \Igl{^) for all e > 0. 

Suppose that Q\i<a]{0) = \Qgl for some > 1, and let \l/ € attain the minimum i"[i,o](e; 0) 
for that /3. Then, 

J [1)0] (e; /3) = £ e)/3 (* e ) = £ e>1 (* £ ) + / S e 2 

>/[i,o](e ) l) + ^ / S e 2 
>i/ M (e) + ^i / 5 e 2 . 

In particular, we conclude that 

U sj->o. 

By the analysis of problem I[ 1)0 i(e; 0) in |ABlj . minimizers have a vortex ball D e (of 
radius 0(e)) with degree deg(?/>+, 83 e ) = 1, deg(?/>l, <9D e ) = 0. Rescaling by e and passing to 
the limit, these converge to a locally minimizing entire solution of ([3]) in M 2 , with degree 
[n + , n_] = [1, 0] at infinity. The convergence being uniform on any compact set, the limit \&* 
has spin 5* = |(|^* | 2 — iV-'+l 2 ) = in M 2 . This contradicts the main result of [ABMj . where 
it is proven that \ip*_ | is bounded away from zero for such solutions (and hence S* > at the 
zero of □ 



We now turn to the case < < 1, where the situation is very different. 
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Lemma 2.3. For all (3 £ (0, 1), 



7T 



g M (/3)<2Q [1)0] (/3)<Q Gi -(l-^)-. 

In particular, for /3 < 1 we may already conclude that the Ginzburg-Landau solution 
is not the minimizer for the problem J^i] (e; /3) , and hence Theorem 11.11 follows from Lem- 
mas GL3] and ETJ It is an interesting open question to determine whether the strict inequality 
Q[i,i]W) < 2( 5[i,o](/3) holds or not. 

Proof. Let u e minimize Igl(c), and set ip + = -^u e , ip- — 77§- Then \1/ = (ip + ,ip^) is 
admissible for /[i,o](e; /3), and 



We;« <£«(*) 



W+ (1+/J) 



16e 2 



1 — \u, 



|2\2 



-GJu f 



(1-/3) 



:i - k 



1 2\ 2 



2" cv " c/ 16e 2 

By blow-up and the result of Brezis, Merle, and Riviere [BMR] we have 

1 



- / (1- \uA 2 ) 2 ^2tt, 



and therefore we conclude that 



Q m {P)<\Q GL -{l-P)\, 



□ 



as claimed. 
Lemma 2.4. For all (3 £ (0, 1), 

QaL<Q[i,i](/3) + (l-/3)limmfi / S 2 dx, 
where S e is the spin associated to any minimizer of Iu\\{e\ 0). 

Proof. Let \E fe minimize I[i,i](e; 0) with < f3 < 1, and set v± = \/2ijj±. Each v± is admissible 
for I GL (e), hence 



I GL (e)<-(G e (v + ) + G e (v_)) 



2e 2 



i,e 1 2 



|0. 



e |2 



(16) 



(1-/3) 



^i,i]( e ;/ ? ) + 



(i-/?) 
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Subtracting n\ In e| from both sides and passing to the limit we obtain the desired conclusion. 

□ 

By putting together Lemma HHH with Lemma |2~41 we obtain interesting information about 
the minimizers of the problem J[i ! i](e; 0): when < /3 < 1, a fixed amount of the core energy 
must come from the spin term. 



Corollary 2.5. If^ e is any minimizer of Jfx i i](e; 0) with < < 1, and S e is its spin, th 



en 



lim inf — / S 2 dx > —. 

e 2 .L 4 



A similar calculation applies also to the fractional degree case: 

Lemma 2.6. Assume < (3 < 1. Then, for any minimizer ty e of iho](e; 0) , its spin S e 
satisfies: 

lim inf — / S 2 dx > —. 

Proof: Let ^f e minimize 7[i t o](e;/3) and v± = y/2ip±. Now only v + is admissible for the 
problem so 

\l GL {e) < \G £ (v + ) 

< ~ [G e (v+) + G t {v.)\ 

< Ee m + Q-Jl [ S 2 



where we have applied the same reasoning in the next-to-last line as in the computation (fl6l) 
above. Subtracting 1 1 In e| from both sides and passing to the limit we have 



2 



Qgl < Q[i,o0) + lim inf 2—^ / S e 2 . 



The conclusion then follows from Lemma 12.31 

3 A different way to measure core energies 



In this section we consider entire solutions ^ in IR 2 , satisfying the integrability condition (jJJ), 
with given degree [n+,n_] at infinity. Measured on all of M. 2 , the energy (defined in (JH])) of 
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such a solution diverges. However, when properly renormalized, there is a well-defined core 
energy, defined as the limit below: 

Lemma 3.1. Let ^ solve ([3]) in M 2 , satisfying (jlj). Then, the following limit exists: 

71 



lim 

R-^oo 



E^;B R )--(n 2 + + n 2 _)\nR . (17) 



For locally minimizing \1/ we expect more. In the case of the single Ginburg-Landau 
equation 

- AU= (1- \U\ 2 )U, (18) 

it is known from Shafrir |Sh] that the only nontrivial solutions which are locally minimizing 
for the energy 



G^to) = JJ±\VU\ 2 + \(\U\ 2 -1) 2 } 



have degree deg(U, oo) = ±1. And for those (unique, by |M2j ) local minimizers, the anal- 
ogous limit ( I17p coincides with the vortex core energy defined via the symmetric Dirichlet 
problem fTTZj) , 

Qgl = lim (G*(u;B R ) -irlnR). 

R-^oo 

We will show that the same is true for \&: 

Proposition 3.2. A nontrivial local minimizer of ([3]) satisfying (j4]) must have degrees n± e 
{0,±1}. 

Proposition 3.3. For any (3 > and for any given degrees n± G {0,±1} ; if ^ is an entire 
solution satisfying (fj]) ; then 



lim 

if,— S>oo 



E^;B R )-^-(\n + \ + \n_\)\nR >Q [n+yT 



2 

If in addition ^ is a local minimizer of energy, then equality holds in the above. 

We begin by proving Lemma 13.11 

Proof of Lemma \3.1\ By the estimates in [ABM] , there exists Ro > for which the solution 
\l/(x) admits a decomposition for \x\ > Ro in the following form: 



ip±(x) = p±{x) exp[ia±(x)], a±(x) = n±9 + x±i x )> 
with x±( x ) <f ) ± (constants) uniformly as \x\ — > oo. 



(19) 



Two-component Ginzburg-Landau 



15 



Without loss of generality we may take 4>± = 0, so X±( x ) ~* uniformly as \x\ — >■ oo. 
Moreover, the estimates in [AB M] imply that for large r, 



P± 



V2 



|Vp±(x)| < 



c 



/ [|Vp ± | 2 + |V X± | 2 ] < oo, 



(20) 
(21) 
(22) 



'\x\>R 

for i? sufficiently large that the decomposition ffT9]) holds. 

First, we observe that for any R > R , by an integration by parts we have: 



Vx± • V0 = 0. 



Thus, 



with 



± 



E^;3 R \B Ro ) -l( n 2 + + n 2 _)ln§- = [ f 

* ^0 JH> R \B Ro 



1 N n ± , 2 IV7 12 . « f 1 A 



I V/^| 2 + ( pi - ^=J ^ + ^|V X± | 2 + 2n ± (p ± - — j • V v± ] 



Using the estimates (T20]) — (T22|) . / is integrable in IR 2 \ H)r ; writing 

E(^;B R )-^-(n 2 + +n 2 _)\nR = E(^;B Ro )-^-(n 2 + + n 2 _)\nR + . 

1 1 Jv>r\Br 

we conclude that the limit R — > oo exists. 

Next, we do a patching argument as in |Shj . 



□ 



Lemma 3.4. Let \I/ fre an entire solution of (0|) satisfying Q). Then, there exists a family 
r E ^(Dr; C 2 ) of functions so that 

^r(x) = ty(x) for \x\ < ~ 7 



i(n + e+0 + ) ; e i(n_e+^_)j ^ on | x | = ^ j or constants ± e M; 



•-R 



|VM/| 2 + (1), / (|^| 2 -1 



(|^| 2 -l) 2 + o(l) 



5 



S 2 + o(l), 
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as R — > oo, where Sr = |(|^_| 2 — \ipn i+ \ 2 ). In particular, 
as R — > oo. 

Proof: We employ the same decomposition ( TT9]) for ip±(x), \x\ > Ro, as in the proof of 
Lemma [3. II Define the cutoff, 

{0 ifr<f, 
32 S£ a , if # < r < iJ, 
1 if r > R. 

We define our modification 

= ^h,+,'0h-J , $h,± = P±(^)exp[ia ± ], 
P±0) = "^| L ( r ) + ( X ~ L ( r ))P±(aO, 5±(ac) = n±9 + (1 - L(r))x±(a?). 
Then, 

|VP ± | 2 - |V P± | 2 = j^- 2 (± - p.)' + (L 2 - 2L)|Vp ± | 2 

By combining the estimates < L(x) < 1, p\ + p 2 < 1, |Vp±| < Cr~ 3 (see [ABMJ) with 
( )22|) . we obtain 



(23) 



|Vp±r- |Vp±r dx->-0 



'H/2<|x|<.R 

as R — > oo. 

Let C =[\n2]-\ Then, 



|Va ± | 2 - |Va ± | a = |V X± | 2 [(1 - Lf - 1] + ~ 2 xl ~ 2L— (V X± ■ 5) 

2C 

-— (l-L)x±(Vx±-r). 



We expand 



pl|Va ± | 2 - pl| V« ± | 2 = p 2 ± (|V« ± | 2 - |Va ± | 2 ) 

+ |Va ± | 2 (p 2 t -p 2 ) 

+ (pl-p|)(|V« ± | 2 -|Va ± | 2 ). (24) 
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Defining A = {R/2 < \x\ < R} and taking each term separately, 
/ p|(|Va ± | 2 -|Va±| 2 ) < / ||Va±| 2 -|Va ± | 2 | 

J A J A 



|Vx±| 2 + ^|Vx±| + ^xi 



< c 



with constant C independent of R. The first term tends to zero by ( 1221) directly. For the 
second, 



A r 



, R dr f 

R/2 r J\x\>R 



1/2 



< 



2vrln2 / |Vx±| 

\x\>R 



1/2 



->o, 



again by (122]) . For the last term, 



R dr 



R/2 1 \x\>R/2 



sup x± — 2tt In 2 sup x± ~~ ^ 0> 



\x\>R/2 



since x± — ^ uniformly as \x\ — > oo. Thus, the first term of (I24j) tends to zero as R — > oo. 
For the second term of ( 124)) . note that by ( 1201 ). 



p| -p±| |Va±| z < 



C n + I—, ,2 

# + IVy±I 2 



-+0, 



again using (1221) for the second piece. The third integral is estimated in a similar way to the 
first one and also vanishes as R — > oo. In conclusion, 



p 2 ± \Va±\ 2 -p±|Va±| 2 



0. 



The estimate ( [20]) also implies 

(|^| 2 -l) 2 + 4S 2 = 0(r- 4 ) and (|^| 2 - l) 2 + 4S 2 = O (r" 4 ) 
Therefore, we have 

f |(|^|2 _ 1} 2 _ ( |^|2 _ 1} 2| + ^2 _ g2 R | ^ Q) 

as — y oo. Putting these results together we obtain 



\E(V R ;B R )-E(*;B R )\ = \E(V R ; A) - A)\ 



as i? — >■ oo, which completes the proof of the lemma. 
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Note that by the same procedure as above, but with the choice 

p±(x) = -j=(l - L R (r)) + L R (r)p±(x), d±(x) = n±9 + L R (r)x±, 

(and L R as in the proof of Lemma 13"^ ) we obtain the opposite patching result, connecting 
a given solution outside a large ball U) R to the symmetric boundary condition on d~B R / 2 : 

Lemma 3.5. Let \I/ be an entire solution of satisfying 0. Then, there exists a family 
^ R G -^(Dr \ Dr/s; C 2 ) of functions so that 





for \x\ = R, 










|V^| 2 = 


-- [ |W| 




Jo R \o R/2 



A , on \x\ = R/2, for constants <f>± G M., 

/ 

[ (i^r-i) 2 =/ (|v&i 2 -i) 2 +o(i), f s R = [ s 2 +o(i) 



JO R \D R/2 

as R — )■ oo, where S R = |(|^_| 2 — |^ i+ | 2 ). In particular, 

E(V R ; B R \ 3 R/2 ) = E(V] B R \ B R/2 ) + o(l) 

as i? oo. 

Proof of Proposition \3. 6 A The proof closely follows that of Theorem 2 in |Shj . Let \1/ be a 
local minimizer. If either \n + \ > 2 or |n_ | > 2, we must have n 2 + + n 2 _ > \n + \ + \n_\, and 
hence Lemma 13.11 implies that for all R sufficiently large, 

£(*;Dp) it, n 9x 7r /t , . . 



By Lemma |33| for i? large we obtain ^ R with constants <f>±, defined in B> R \H R /2- Denote 

by 

G,(U;Q):=G t=1 (U;Q)= [ !h W| 2 + \(\U\ 2 - l) 2 



2' 1 4 

the Ginzburg-Landau energy for U G Hl oc (Q;C). Taking Q = D fl/2 , let C/± minimize the 

i/2 



Ginzburg-Landau energy with boundary condition U R \gn R/2 = e l [ n±9+9i± ]. By the results 



of Brezis, Bethuel, & Helein [BBH] . 

©«/ 2 ) = 7r|n±| \n{R/2) + 0(1). 
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Now let 

t*ii(a;), inDfl\D fl/2 
Since \P is a local minimizer, we have 

E{V) B R ) < E($ R ; Br) 

= l - [G*(U+; B R/2 ) + G*(U r] B r/2 )} + E^ R ; B R \ B R/2 ) 
= £ (Kl + \n-\) ln(R/2) + E(V; B R \ B R/2 ) + 0(1). 



From Lemma [3. 1[ it follows that 

Efif] Br \ Br/ 2 ) = |(4 + n 2 _)ln2 + o(l), 
as R — > oo, and hence 

E(V;B R ) < |(|7i + | + |n_|)lni2 + 0(l). 
This contradicts (1251). □ 



Proof of Proposition Iff. 31 Let \l/ be an entire solution to (|3]) satisfying (j!J) with degrees n± G 
{0, 1} at infinity. (The degree -1 case may be obtained by complex conjugation.) Multiplying 
each component ip± by a complex constant of modulus one if necessary, we may assume that 
if)± — > -^e in±9 as \x\ —> oo (that is, (f>± = 0.) For large R, let ^ R be as in Lemma l3^| so 
that ^R\ao R = -^(e in+e ,e in - e ). Note that by scaling, 

J [n+ , n _]{l,R\P) = mm^E(*;B R ) : tf| SDfl = ^(e m+e , e in - e ) J = J [n+>fl _] . (26) 

Using (TP2"]) we conclude that 

E(V;B R ) = E($ R ,B R ) - o{l) 

71 

> J[n+,n-}(l,R] p) -o(l) = ~(n+ + n_)ln J R + Q [ri+in _] - o(l), 
which proves the first assertion in Proposition 13.31 



Now assume in addition that $ is a local minimizer. We now employ Lemma 13.51 with 

= (i>R,+,i>R,-) , $r,± = p±(x) exp[i&±], 



P±(x) = -^=(1 - L R (r)) + L R (r)p±(x), &±(x) = n±6 + L(r)x±{x) 



(27) 
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so that 

^n\do R/2 = -j=(e tn+e ,e m ~ e ), ^B.\ao R = ^\ao R , 

and for all R > R®, 

E($ R ; B R \B R/2 ) = E^ ] 3 R \B R/2 ) + o(l) = ^(n + + nJ) In 2 + o(l). 

Next, we let Vr(x) be a minimizer for the problem I[ n+)1l _](l, -R/2; /3) with symmetric data. 
By scaling, 

E(V R ;B R/2 ) = J [n+ , nA {l,R/2-p) = J [n+ , n _ ] (2/ J R;/3) = |(n + +n_)ln f ^ +Q [n+ ,„_] + o(l). 
We now define a trial function in 3 R , 

&r{x) ■ 

Since § R = \I/ on <9Dr, and \I/ is a local minimizer, we have: 



* fl (ar), if < |ar| < i? 
if |x| < i2/2 



E(*;B R ) < E($ R ;B R ) 



E{V R ; 3 R/2 ) + D« \ D fl/2 ) 

7T 

-(ra + + n„)ln J R + Q [n+in _] + o(l, 



□ 



4 Locally minimizing solutions 

We are ready to prove Theorem II . 31 stated in the Introduction, concerning locally minimizing 
solutions of the system (j3J). 

Proof of Theorem \1. Si. First, if \E r * has the form 

^* = -j= (u(x - a)e l<t,+ ,u(x - a)e^-) , 



then it is a local minimizer. By Theorem the minimizer for the boundary- value problem 
Jn ) i](e; (5) coincides with (v^2 times) the Ginzburg-Landau minimizer in each component. By 
blow-up around the vortex we converge to a solution to the system in M 2 which is a local 
minimizer. Since each component of the minimizer in Di is a solution of the Ginzburg- 
Landau equation, the solution we obtain is of the form above, and hence such solutions are 
local minimizers. 
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Now, assume \&* is any local minimizer with degree deg(/0±, oo) = 1. We claim that its 
spin 5* = |(|V'-| 2 - IV'+I 2 ) = 0. Suppose not, and let 

a = I Sldx> 0. 



Choose Rq large enough so that S^dx > a/2. 

We use Lemma 13.41 with n± = 1 to bound the energies E(^*;B R ) from below. Let 
= [ipR,+,i>R,-] be as in Lemma 13^41 u R) ± = \/2?p R ,±, and apply the identity (|T5l) : 

S(**;Dfl) = S(* fl ;Dfl) + o(l) 

4(aK + ;D fl ) + G»K-;D fl )) + (^-l) / ^ + o(l) 

>Igl + 

= 7T In i? + Q GL + {p - 1)| + o(l), (28) 
since li^i = e* e on <9Dr, and by scaling we have 

J GL (1,R) : = min{a(n;D fi ) : u e = I G l(1/R). 

Now let be as in Lemma 13.51 and define 



Vr(x), ifR/2<\x\<R, 
75 K/ 2 e^+,u R/2 e^-) , if |x| < R/2, 



where ur/ 2 is the minimizer of the Ginzburg-Landau functional in B R / 2 with symmetric 

G*(u R/2 ;3 R/2 ) = mm{G*{v;B R/2 ) : v G H l ei e{B R/2 )} = I GL (2/R). 
By the estimates of Lemma 



E{$ R - B R ) = E$ R ; B R \ B R/2 ) + E($ R ; B R/2 ) 

= E(V*; B R \ Br/2) + Igl{2/R) + o(l) 

^ ln f "7777T J + o(l) J + (tt ln(i?/2) + Q GL ) + o(l) 



= 7rlni? + Q GL + o(l), 
as i? — > oo. Hence, comparing with ( |28|) we have 



fl ; © fl ) < £(^*; D fl ) - (/? - 1)| + o(l), 
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and for a > we have that \&* cannot be a local minimizer for f3 > 1. Therefore the claim 
is established, 5** = for any local minimizer. 

We then conclude that for any local minimizer, |-?/>_| = So ip± have their zeros 

coincident, and each solves 

-A^ = (l-2|^ ± | 2 )^ ± 

in R 2 with degree one at infinity. By Mironescu |M2] each is (a/2 times) an equivariant entire 
solution of the Ginzburg-Landau equation, and therefore \l/ must have the form (jH]). This 
completes the proof of (i) of Theorem 11.31 

Next we turn to the case < (3 < 1. We first show that if \P is a local minimizing entire 
solution in M 2 with n± — 1, then 

f S 2 dx>-, 

with S= |(|^_| 2 - |^+| 2 ). Let 

ffR '■= / S 2 dx, a = I S 2 dx. 

First, we bound E(ty;"Bn) from below. Apply Lemma 13.41 to obtain for each R > a 
corresponding and set u± = ip R ±. By Lemma [33] and (1151) . 



= £(*fl;B fl )-o(l) 

= i(a(« + ;Dfl) + a(w-;Dfl))-(l-/3) / S 2 -o(l) 

>Igl(1/R)-(1-P)(Tr-o(1) 

= Tx\nR + Q GL -(l- f3)a-o(l). (29) 

For the upper bound we use the local minimality of First, as in |BBH] (see [ABlJ) 
for any pair of points ai,<Z2 G Di, p > fixed, and all small e > there exists $ with 
®Ui = Ti^VI so that 

E e (k) = 7T In Qj + W (a l5 a 2 ; [1, 0], [0, 1]) + I [1)0] (e; p) + I [0>1] (e; p) + O(p) 

-lnQ)-| [ln(l - |ai| 2 ) + ln(l - |a 2 | 2 )] 
+ ^[i,o](e,p; /3) + J[o,i](t,p] P) + 0(p), 

where W is the renormalized energy associated to the Dirichlet problem (see |ABlj .) For 
any r\ > 0, fix a\ = — a 2 ^ close enough to the origin and p > such that 

£ e ($ e ) < vrln + J [li0] (e,p;/3) + J [0 ,i](e,p;/3)+r/. 
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Using scaling and Lemma [2.31 



E e (6 e ) < 7r In + J [1)0] (e, p; f3) + J m (e,p;f3)+n 



+ 77 + 0(1; 



< tt| lne| + Qgl - (1 - /3)^ + 77 + o(l), (30) 

where o(l) is with respect to e — > 0. 

We now construct a function with the same boundary values as \I/ on <9B#. Applying 
Lemma[331 we set $ R = V R in B R \B R/2 , with = ^=( e ie+ *+ ,e ie+ *~) } <j>± constants. 

Inside D fi / 2 we define e = 2/i? and $ R = $ 2 /fl(2x/i2), with $ e as constructed above. By 
rescaling and Lemma [3.51 we have (as R — > 00), 

E($ R ; B R ) = E($ R ; B R \ B R/2 ) + E($ R ; B r/2 ) + o(l) 

= 7rln (^J + #a($2/Ji)+°( 1 ) 

< 7T In/? + Q GL - (1 - $)j + 77 + o(l). 

Since \1/ is a local minimizer, we obtain the upper bound, 

E{V; Br) < vr In R + Q GL - (1 - 0)^ + v + o(l), (31) 

for any 7/ > 0, in the limit R — > 00. 

Putting together ([29]) and f[3Tj) . since /3 < 1 we have 

7T 77 
(7 > - 



4 1-/3 

Since 77 > is arbitrary we obtain the estimate (J7J) in Theorem 11.31 As a corollary, we 
observe that the equivariant, Ginzburg-Landau-like solutions are not local minimizers for 

< /3 < 1. 

To complete the proof of Theorem II .31 we require the following uniqueness result for entire 
solutions with a common zero in each component: 

Proposition 4.1. Suppose ^> is an entire solution of (j3]) ; satisfying (j4j), with [n + ,n_] = 
[1,1]. If in addition tp + (xo) = ip-(xo) = for some xq G M. 2 , then there exist constants 
(f>± e M with 

I 1 

ip+(x) = —=u{x - x )e^+, = - x )e^-, (32) 

v 2 v 2 

where u = f{r)e ld is the equivariant, degree one entire solution to the Ginzburg-Landau 
equation in IR 2 . 
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Remark 4.2. The proof of Proposition 14.11 is based on the method of Mironescu |M2j . and 
is deferred to the end of this section. The existence or non-existence of a locally minimizing 
entire solution for fl3]) with degrees n + = 1 = n_ is an interesting open question. If such 
a local minimizer were to exist, Proposition 14.11 implies that its vortices must be spatially 
separated. 

We may now complete the proof of Theorem II .31 Since a > |, the local minimizer \l/ 
is not of the form ( 1321) . and hence by Proposition 14.11 we conclude that |\I/(x)| 2 ^ on I 2 . 
By (l20il there exists R > so that |\l/(a;)| 2 > \ for |x| > R, while on the compact set Dr, 
|\I/(:e)| 2 is continuous and non-vanishing, thus bounded away from zero. This completes the 
proof of Theorem 11.31 □ 

Remark 4.3. When (3 = 1 the equations decouple completely, and each component satisfies 
the Ginzburg-Landau equation in M 2 . In this special case there is a huge degeneracy, and 
the solution space with degree [n +, nJ\ = [1,1] is completely described by: 



where u = f{r)e %9 is the unique equivariant solution with degree one, a,b el 2 are arbitrary 
points in the plane, and <j)± are arbitrary constants. Since the equivariant solution is a local 
minimizer for Ginzburg-Landau and because the energy decouples when (3 = 1, 



it is clear that each of these solutions is a local minimizer. The question then arises whether 
there exist local minimizers with (3 < 1 which bifurcate from this degenerate family as /3 — > 1. 
Alas, if we calculate the L 2 norm of the spin S a ,b,<t>+,^- associated to C/ a> 6 j( ^ +j ^_, we have 



and the value approaches 7r/4as|a — b\ — > oo. In particular, by Theorem 11.31 no family ^ p 
of local minimizers for (3 1~ can converge to a member of the family U a ^ + ^_, and so 
there is no bifurcation. 



U(x) = U a ,b,</> + </> 



-={u{x- a)e**+ ,u{x- b)e*- ) , 





dx < — , 
4 



Next, we complete the proof of Theorem 11.41 stated in the Introduction: 
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Proof of Theorem LJ., Let e Hg(Q\ C 2 ) be minimizers of E £nj p with < /3 < 1. We argue 
by contradiction, and assume that there is no such constant c > for which |^| > c in Q. In 
that case, there exists a sequence e n — > and points p n G ft for which ^/ en (p n ) — > 0. By the 
methods of [BBH] extended to the energy E £ (see [AB2] .) |^ n | is bounded away from zero in 
a neighborhood of the boundary dfl, so in particular, we can assume dist(p n , dQ) ^> e n for e n 
sufficiently small. Blowing up around p n at scale e n , by standard estimates a subsequence of 
the rescaled minimizers converge in Cf oc to a locally minimizing solution \I/ to ([3]), satisfying 
(HI), with \I> (0) = 0. By Proposition 13.21 the degree pair deg(\I/, oo) = [ra+,n_] has one of 
the forms [0,0], [±1,0], [0,±1], [±1,±1]. We claim that each is impossible. Indeed, from 
the only [0,0] solution is constant (of modulus l/y/2) in each component. In [A.BMJ 
it is proven that the [0, ±1] and [±1,0] locally minimizing solutions never vanish for any 
(3 > 0. And Theorem 11.31 (ii) asserts the same conclusion for the [±1,±1] local minimizer 
when < /3 < 1. In conclusion, there must exist c > for which |\EU > c > as claimed. □ 



We conclude the section with the deferred proof of Proposition 14.11 



Proof of Proposition 4jJ_ ■ We follow the uniqueness proof of Mironescu [M2] . Denote by 
P—^f- Without loss of generality, assume xq = 0, and form the quotients 

^) = ^6C»(R a \{0}). 
We may then derive the system of equations satisfied by w±, 



-Aw + - 2 
-Aw_ - 2 



p' <9w + 
p dr 
p' dw_ 



1 



— w + = p 2 [2 - \w + \ 2 - |w_| 2 + /3(|w_| 2 - \w+\ 



1 



—w- = P 2 [2 - \w+\ 



P(\w-\ 



w. 



(33) 



p dr r' 

As in the derivation of the Pohozaev identity, we multiply each equation by x ■ Vw±, and 
integrate over the domain Djj \ O# . The resulting identity has the form 



with 



F(R)-F(R ) = G(R,R 



F(R) 



Rp 2 



[(2-\w + \ 2 -\w_\ 2 ) 2 + P(\w-\ 2 -\w + \ 2 ) 2 } 



E 

± 



^k±T + y [\d r w±\ 2 - \d T w±\ 2 ] 
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and 

'2rp' 



G(R, Rt 



(\d r w + \ 2 + \d r w^\ 2 ) 

[(2-K| 2 -| w _| 2 ) 2 + /3(k-l 2 -KI 2 ) 2 ] 



p 

2p 2 + 2rpp' 



4 

Using the estimates (|20p - (|22p on the solutions ift± for large R, the corresponding facts for the 
Ginzburg-Landau profile p, we obtain F(R) — F(Rq) — > as both R — > oo and Rq — > 0, and 
thus conclude that G(oo, 0) = 0. Since p > and p' > 0, we may conclude that \w+\ = \w~\ = 
1 and d r w± = in all of R 2 , and hence w± = w±{6) = e %x± ^\ that is ip± = pw± = pe lx± ( e >. 
By the estimates in |ABMj . if)± — -^e 1 ^^ — > uniformly as |x| — >■ oo, with constants 
0-1- G JR.. In particular, x± = 4>± are constant, and ip± = pe^ e+x± ^ = -±= J(r)e^ e+9i± - ) as 
claimed. □ 

5 Bifurcation of symmetric vortices 

In this section we study the stability and bifurcation of the equivariant solutions of the 
Dirichlet problem in the unit disk Di. For convenience, we replace the usual parameter e by 
A = e~ 2 in both E e and the Ginzburg-Landau energy G e , and (with abuse of notation) write 

Ex(?) = I {^|V^| 2 + ^(|v&| 2 -l) 2 + M ( |^ + | 2 _|^| 2 ) 2 |^ 



G\(u) 



I| Vm | 2 + ^(| m | 



We consider critical points * G H := G H l (Bi,C 2 ) : $?\ mi = ^{e ie , e i6 )} , which solve 
the Dirichlet problem, 

-A^ + =A(1-|^| 2 )^+ + A/3(|^-| 2 -|^+| 2 )^+ inBi 
-Aip_ = A(l - |^| 2 )^_ - \(3(\ip-\ 2 - \ip+\ 2 )ip- inDi 



i>± = -^-=e iB on 9Di 

v2 



(34) 



in the unit disk Di, for fixed (3, as A ranges in the half-line A G (0, oo). We will show that if 
< (5 < 1, the symmetric vortex solution 

where u\ = f\(r)e z6 is the degree-one equivariant solution of the classical Ginzburg-Landau 
model in the unit disk, is stable (a strict minimizer of energy) for A = e~ 2 small, but loses 
stability at some A^ = e^ 2 > 0. 

First, we note that for A small enough, there are no other solutions to 
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Proposition 5.1. There exists A^ > so that for every A < A^ the unique solution to ( 1341) 
is the symmetric solution U\. 

Proof. First, define the convex set B = G H : \^(x)\ < 1 in O x }. By Lemma 4.2 of [AB2j . 
any solution of (|34|) lies in £?. The second variation of the energy around \l/ = (?/> + , f_) 6H 
in direction $ = G /^(©i; C 2 ) is: 

= / {|V$| 2 + A(|^| 2 -1)|$| 2 + 2A(^,$) 2 

+/3A [|^ + | 2 - \^-\ 2 } [\<p+\ 3 - l^-l 2 ] + 2/3A [<^ + ,</> + > - (^-,^)] 2 } • 
For any \I/ G £>, we have 

^($)[$]> / {|V$| 2 -C(/3)A|$| 2 }rfx, 

with constant C(/3) > independent of A, $. By choosing A^ sufficiently small that C((3)\ 
is smaller than the first Dirichlet eigenvalue of the Laplacian in ©! we may conclude that 
E'K^f) is a strictly positive definite quadratic form on Hq(JD>i; C), for any $ E B. Thus, -Ea 
is strictly convex on B, and hence it has a unique critical point. □ 

We also observe that the symmetric vortex solution U\ is the unique solution to ( 1341) for 
which both components vanish at the origin: 

Proposition 5.2. Suppose ^ is a solution to ( 1341) with ip+{0) = = ip-(0). Then \I/ = U\. 

Proof. The proof is exactly as for Proposition 14.11 as before, let w± = Jrfiy> with p(r) = 
^|«a|, so w± = e td on dTS>i. In particular, \w±\ 2 = \d T w±\ 2 = 1 on <9Bi, so 

F(l) = -\ I \d r w+\ 2 + |<9 r w_| 2 < 0. 

Arguing as above, we may conclude that (7(1,0) = 0, and hence \w±\ = 1 and d r w± = in 
Di, with w± = e ie on dH>i. Thus, w±(x) = e %e on D 1; that is ip± = p(r)e ie = -j^u\. □ 

It is important to recognize the role played by two groups of symmetry acting on the 
problem. First, the problem is invariant under the action of the group § , with the following 
representation: for ip G H ie (Di; C), writing the independent variable in complex form z = 
x + iy, we represent £ G S 1 C C via 



(R^)(z)=^z). 
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By abuse of notation, for \1/ = G H we define 

= (R^+ , R^-) . 

Our problem is also invariant with respect to the involution, 

TV — (-V-(-x) , —v + (—x)) , V = (v+(x),v-(x)). (35) 

The equivariant solution U x is fixed by both groups, R^U X = U x = TU\, for all £ G S 1 . 

With the goal of studying bifurcations from the symmetric solutions, our first task is to 
study the spectrum of the linearization of the energy around U x , 

&x(U x )[*]= [ [|V$| 2 + A(K| 2 -l)|$| 2 + A( MA ,^ + + ^-} 2 + A/3K,v9 + -^) 2 ], 

for $ = (</?+,</?_) G Hq(Di,C 2 ). The quadratic form E"(U X ) is associated to the linearized 
operator, 

-Ay? + + A(|m a | 2 - l)v+ + A(m a , </?+ + v?_)m a + /?A(m a , y? + - y?_)w A 
-Ay?_ + A(|m a | 2 - IV- + A(m a , y?+ + V?-) m a - /?A(m a , y?+ - y?_)w A _ 

For each A > 0, L x defines a self-adjoint operator acting on its domain H 2 n ^(Di; C 2 ) C 
L 2 (Di;C 2 ) (with real scalar product). Moreover, by elliptic regularity theory, for all A > 
is has compact resolvent, and thus discrete spectrum, consisting of eigenvalues, 

<t(La) = WA) < MA) < MA) <•••}, 

repeated according to their (finite) multiplicities, and ordered by the min-max principle for 
each fixed A > 0. . We seek the critical value A = \p at which the radial solution loses 
stability with increasing A. The group invariance under R^, £ G S 1 , and T are inherited by 
the linearized operator, and each of the eigenspaces is invariant under these two symmetries. 

We will uncover the spectral properties of the linearized operator L x via a sequence of 
reductions, and prove analytic dependence on A of eigenvalues and eigenfunctions provided 
they are simple modulo the action of the group S 1 . 

Reduction of L\ 

First, we reduce to a scalar problem: let 

C X (p : = -Ay? + A (|M 2 - 1)<P + 2/3A (u x , (p)u x , 
as an operator acting on L 2 (Di; C), with associated quadratic form 

Q x {w):= [ [\Vw\ 2 + \(\u x \ 2 -l)\w\ 2 + 2(3\(u x ,w) 2 }. 







L-ct>_ 
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We also define the linearization of the classical Ginzburg-Landau energy, 
Gl(u x M= [ OVyf + A(M 2 -1)M 2 + 2A( MA ^> 2 ]. 

and 

Cfip := -Aip + A(|w A | 2 - l)ip + 2A(w A , ip)u x . 
By |Mlj . is a positive definite operator for all A > 0. 

Lemma 5.3. (1) For any > 1, L x is positive definite for all A > 0. 

(2) For any (3 > 0, \x is an eigenvalue of L x with eigenfunction $ if and only if: 

(i) either \i 6 cr(Cf); 

(ii) or // G ct(£a) and $ = (if, —(f), where <p is an eigenfunction for C\. 
Proof. A simple calculation shows that 

E'l(U x m = Gl(u x )[cp + ] + Gl(u x )[^} + ^-l)X [ (u x ,v + - V _) 2 . (36) 

By [M1J. G'{(u x ) is positive definite, and thus (36) implies that L x is positive definite for 
P > 1, proving (1). 

Now let (3 > be arbitrary, and assume $ = (<p+, <p~) solves L A $ = We then have 

L+$ = C%<p+ + (/3- 1)A (u x , u x 
L A $ = Cfip- -{p- 1)A (u x , (p+-(p-) u x . 

Adding these two identities together, 

Cf(ip + + </?-) = L$$ + = /j(ip + + tp_). 

Thus, if /i ^ a(C^), we must have — (p- = (p + =: y?, and moreover L^Q = £\(p = fi(p. 

This proves the only if part. The "if" part is obvious. □ 

Denote the eigenvalues (repeated by multiplicity) of £ A by 

*(£x) = UkM < <•••}• 

The next step is to decompose w G Hq(D; C) in its Fourier modes in 9: 
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Using Parseval's identity, we have 

/ {u x ,w} 2 dx = - flir) \b n+ i + b 1 ^ n \ 2 rdr. 

Consequently, the operator C\ can be identified to a direct sum in Fourier modes, 

oo 

£ AW ^04 n) (6 n+1 A 



>\-n) 



(37) 



n=0 



where the operators C^/ are associated to the quadratic forms 



Jo 



{n + If 



:i-n) 



\K+l\ 2 + \b'l- n \ 2 H 5 l^n+ll 2 5 l^l-n| 2 



+ A(/ 2 - 1) (|& n+1 | 2 + | 0l _ n | 2 ) + XI3fl\b n+l + 6i_ n | 



for n^O, and 



r dr. 



If we let 



X = {/ : (0, 1] C; /(r)/r G L 2 , /' G L 2 , /(l) = 0}, 
y = {/:(0,l]^C;/GL 2 ,/(l) = 0}, 

then: acts on X, Q ( p acts on X © F and, for n ^ 0, 1, acts on X © X. 

The spectrum of this direct sum is given by the union of the eigenvalues of the operators 
. As we show below, c 
Define (in Fourier Space) 



£^ . As we show below, only one of these operators may contribute an eigenvalue near zero 



£\w := (J)£i n) (6 n+ iA_ n ), 

where w = ^2 n ^ 2 b n (r)e md , and so £\w = C^\b 2 ,b ) © C\w. Let Q\ denote the quadratic 
form associated to C\. 

Lemma 5.4. There exists a function rao(A) > such that inf o~(£\) > mo(A) > for every 
A > and, in addition, inf mo(A) > /or eac/i compact interval I C (0, oo). 

Proof. This follows from the reductions described on page 337 of |Mlj . Let w G Hq(Dx; C) 
with associated coefficients {o«(r)} n ^o,2, such that Xln^o 2 IIMIi 2 ((o lWrfr) = ^ e define 
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°i := ^(Zln^o,2 l fo n| 2 ) 1/2 - Then, ai(r) is purely imaginary, ||ai|| L 2 ((01) . rdr) = 1, and 
Q x (w) > tt E f (|&' n (r)| 2 + ^|6„(r)| 2 + A(/ 2 (r) - l)|& n (r)| 2 ) r dr 



n^0,2 
-I 



> ^ I \ KM! 2 + + A(/ A 2 (r) - l)h(r)| 2 \rdr 



= ^Qi 0) («i>«i)- 

Now set TOq(A) := inf iQ^^ax, ai), where the infimum is taken over all d\ G Hq((0, 1)) with 
|| a i||L 2 ((o,i);»*rfr) = 1- in [Ml] it is proven that to (A) > for all A > 0. By the min-max 
principle, inf cr(£ A ) > m (A) > 0. 

It remains to show that mo (A) can be bounded away from zero when A is bounded 
away from zero. To do that, we make a transformation as in [Ml] in order to consider a 
fixed system of equations on an increasing family of disks Dr, with R = a/A. We define 
Ur(x) := Ux(x/y/\) = F(r,R)e td , and di(r) = a x (r/v 7 !), for r G [0,R]. In this way, 



m (A) 

^ H a ilL2 ((0 , fl) . rdr) =i 



m (i?) := n _ m inf tt / <j | | 2 + — a 2 (r) + (F(r, R) 2 {r) - l)|ai| 2 [> r dr. 



We now show to (-R) is decreasing in > 0. First, the radial profiles /r(t*) are pointwise 
decreasing in i?: whenever R < R', we have F(r, R') < F(r, R) for all r G (0, R). (Indeed, on 
[0, R], F(-, R') is a subsolution for the equation satisfied by F(-, R).) Using this fact and the 
inclusion Hq(3r) C Hq(Br'), we find that rho(R') < rho(R). In case the two are equal, this 
would imply that the minimizers of in (which are nonnegative,) vanish identically in 
3ri \ 3r. This is impossible, by the maximum principle. Thus m (-R) is strictly decreasing. 
Finally, given a fixed interval I = [A , Ai], whenever A G I, 

Wo (A) = Am (v / A) > A to (\/Ai) := Cj > 0. 

□ 

Combining the results of Lemmas 15.31 and 15.41 we may conclude: 

Corollary 5.5. If /2 k (X) G cr(L\) and //fc(A) < for some A > 0, then /^(A) G a(C^). 

In particular, if yUi(A) is to cross zero it must be because of the ground state eigenvalue 
/4 (A) of Z^. We perform one final reduction of the operator £ A . We define a quadratic 
form Q A on real-valued radial functions (ao, 02) by 



Qa (00,02) := tt / 
./o 



i'o) 2 + (4) 2 + ^«2 + A(/ A 2 - 1) (o 2 + a 2 ) + (3\fl (a - a 2 ) 



r dr. 
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The self-adjoint operator associated to Q x is 



«0 
«2 



-eft - -a' + X(fl - l)ao + 0\fl(a o - a 2 ) 
1 4 

-d2 - -a' 2 + — a 2 + A(/x - l)a 2 - WaK - a 2 ) 



(38) 



Lemma 5.6. /i G I is on eigenvalue of C x over L 2 (([0, 1]; rdr); C 2 ) if and only if it is an 
eigenvalue of Ai\ over L 2 (([0, 1]; r dr); M. 2 ). Moreover, if fi is a simple eigenvalue of Ai\ 
with eigenspace spanned by (00,02), then 



ker(£$ i 1) 



fil) = {t(£a ,-£a 2 ): £ e S\ teR}. 



Proof. Let /i G a(C^) with (complex- valued) eigenvector (b ,b 2 ), that is: 



C 



(i) 



-60 



6 + A(/ A 2 -l)6o + W 2 (6 -& 2 ) 



-&2 " ~&' 2 + 4&2 + A(/ 2 - 1)62 + /3A/ 2 (6 2 - 60) 



(39) 



We observe that a = Im6 , ^2 = Imfe 2 will be eigenvectors of M.\ with the same eigenvalue 
/j,. On the other hand, it is clear that if (a , a 2 ) are (real-valued) eigenvectors of then 
(b ,b 2 ) = (ia ,ia 2 ) will be eigenvectors of C x ^ with the same eigenvalue. Thus, cr(C^) = 
o-(M x ). 

Finally, suppose /z is a simple eigenvalue of J\4\ with eigenspace spanned by (a ,a 2 ). If 



(&o, b 2 ) is an eigenfunction of C x , then (by the observation above) (Im6 , Im& 2 ) 
for g G R. Similarly, (Reb Q , —Reb 2 ) is an eigenfunction of A4\, and so (Rebo, —Reb 2 ) 



-q{a ,a 2 ) 



p(ao, a 2 ) for p G M. Setting t = ^/ p 2 + q 2 and £ 
t (£a>o , — £ a 2 ) , as claimed. 



p + iq 
~~t 



G S , we have (b ,b 2 ) 



□ 



We conclude this part with the following essential fact about the ground state eigenvalue 
of My. 

Lemma 5.7. The ground state eigenvalue n\ (A) 0/ .Ma *5 simple. It is generated by 
(a (r; A), a 2 (r; A)) with < a 2 (r;A) < a (r;A), a (r; A), a 2 (r; A) > m (0,1), a (0;A) > 0, 
and a 2 (r; A) = 0(r 2 ) /or r — > 0, for all A > 0. 

Proof. This fact follows as in |M1] : we claim that, up to a change of sign, < a 2 (r) < ao(r) 
holds for all r. Indeed, if not we define d 2 (r) = min{|ao(r)|, |a 2 (r)|}, and similarly do(r) = 
max{|a ( r )|? l a 2(^)|}- Replacing a , a 2 respectively by So, d 2 does not change the quantity 

/ [| a o( r )| 2 + \a 2 (r)\ 2 ]r dr, and the first, second, and fourth terms in Q x ^ are unchanged. 
Jo 
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However, the third and last terms are reduced, contradicting the minimality of the Rayleigh 
quotient at (00,02). By a standard argument, nonnegativity of the eigenfunctions and of 
ao — 02 implies simplicity of the eigenvalue. Each function is strictly positive in (0, 1) by the 
strong maximum principle (or the uniqueness theorem for ordinary differential equations). 

The behavior of do, a 2 at r = follows from the ordinary differential equations (see (16) 
of |Mlj ) satisfied by a , a 2 in (0, 1). □ 

Analyticity 

In this part, we prove that the radial profile f\ is analytic in both A and r, and conclude 
real analytic dependence on A of the simple eigenvalues of L\. 

Proposition 5.8. f\(r) is real- analytic in r G [0, 1] and A > 0. 

As in the proof of Lemma 15.41 we rescale our problem to study a fixed equation in a 
variable domain. We define 

F(r,R) = f x (r/VX), 
for r G (0, R), with R = Then, F solves 

- F" - -F' + i-F = (1 - F 2 )F, F(0) = 0, F(R) = 1. (40) 

While the existence of such an F(r) may easily be done by minimization of its energy 
functional, to obtain the desired properties of F it will be necessary to relate the solution 
F(r, R) of the boundary value problem to solutions <j>(r, b) of a Cauchy-type problem, 

- <t>rr ~ -<j>r + k<P = (1 - 0(0, 6) = 0, 0'(O, 6) = 6, (41) 

where the value of b is chosen (by "shooting") to achieve the boundary condition <p(R, b) = 1. 
Indeed, it has been shown (see |CEQ]) that F(r, R) = (f)(r, bji) for a unique value of b = bn > 
0. We note that the equation being singular at r = 0, this is not a regular initial- value 
problem, and thus the existence and analyticity of the solution do not follow directly from 
the Picard existence theorem. (See Theorems 8.1 and 8.3 of |CL] .) 

In the remainder of this section it will be convenient to extend <p(r, b) to r G C, with 
complex parameter b G C. The following equivalence follows easily from the variation of 
parameters formula, and may be established by direct calculation. 

Lemma 5.9. Let p(r) be continuous for \r\ < r , r > 0. If f is continuous on \r\ < r and 
solves 

/( r ) = 6r+- / (i-L\ S g{s)ds, (42) 
^ Jo x - r s/ 
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then f G C 2 ({r G C : < |r| < tq}) and solves 

-frr--fr + \f = 9(r), forre(0,r ), /(O, 6) = 0, f'(0, b) = b. (43) 

Conversely, if f G C^D^) U C 2 (O ro \ {0}) is a solution of ([4"3~]1 . i/ien a/so so/wes fl42l) . 
Moreover, if g is analytic in D ro , then f will be analytic for r G D ro and b G C. 

The integral in (l4"2{) is a complex path integral, and the path is the straight-line segment 
joining to r in C. Note that in the analytic case, the singularity of the solution of fj4*3|) is 
removable, and thus / is analytic in the entire disk D ro . 

Using Lemma 15.91 we may obtain an analytic solution to the Cauchy-like problem ( l4Tj) 
by Picard iteration of the integral equation, 

0( r , h )= hr+ ( s - - L) 8 (1 _ 0( s , h f)<P{s, b) ds. (44) 
Let b G C be fixed, and define a rectangle in C x C, 

1 1 



1Z := {(r, 6) : \r\ < p, |6 — 6q| < p} with p := min 



2' 2|6 n 



Lemma 5.10. For any fixed bo G C, there exists a unique solution to fT4T]) . which is analytic 
for (r, b) E 71. 

Proof. Let h(r,b) = <p(r,b)/r. Then, / solves (|44p if and only if h solves the fixed-point 
equation, 

h{r, b) = T[h]{r, b) := b + i jf ^ - 1^ s(l - s 2 /i(s, 6) 2 )/i(s, 6) ds. 

Define a class of continuous functions, 

M := {/i G C(D p ; C) : \h(r) - b \ < 1 Vr G © p } . 
We observe that if /i G M, then for any r G B p , 

|Wi(r)| < |r(/i - Ml + \rbo\ < ^ + ^ = 1- (45) 
Thus, |1 - s 2 h 2 (s)\ < 2 for < |s| < |r| < p, and 

|T[/i](r) - fool < \b - bo\ + \ /' |1 - s 2 /i 2 (s)| d|s| < J + |r| < 1, 

2 Jo 2 
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and hence T : M -»■ M. If hi, h 2 G Af, then 

|T[/i 2 ](r) - T[/n](r)| < J [l + |s| 2 (M 2 + |/ii/i 2 | + |/i 2 | 2 )] |s| M«) - hi(s)\ d\s\ 

^ Jo 

r\r\ 

< 2\\h 2 - 1 00 / k| d|s| 

= |r| 2 \\h 2 - /ii| I do < ^||/i2 - /ii|U, 

for any r G D p , where we have used (]4"5]) to estimate the integrand. Thus, T is a contraction 
on Af , and there exists a unique continuous solution of the fixed point equation for r G D p . 
Since the solution may be characterized as the uniform limit of the iterates h n := T[h n -i], 
ho = b, and by induction each h n is analytic for (r, b) G 1Z, the solution h so obtained is 
analytic for (r,b) G Setting 4>(r,b) = rh(r,b), we obtain a unique analytic solution to 
f !44p . which by Lemma [5.91 provides a unique analytic solution of (141 j) for (r, 6) e 7£. □ 

For r ^ 0, the initial value problem for the differential equation ( 14 ip is regular, and the 
existence, uniqueness, and analytic dependence of solutions in the complex plane follow from 
Theorem 8.3 of [CL] . Using this observation, we may extend the solution 0(r, b) analytically 
along the real axis in r to obtain a maximally defined analytic solution. Indeed, for any 
b, \b — b \ < p, define Rb to be the supremum of all real values R > for which there 
exists an analytic solution 0(r, b) of ( 141"]) in a C-neighborhood of the real interval [0, R]. By 
Lemma [5. 10] we know that Rb > p. 

Proposition 5.11. For any b G M there exists a unique solution </>(r, b) to (|4ip forr G [0, Rb), 
which is real analytic in (r, b). Either Rb = 00, or Rb < 00 and lim^^- |0(r, b)\ = 00. 

Proof. The proof is a standard extension argument from the theory of ODE (see |CL] .) By 
the definition of Rb, for each real value s G [0, Rb) there exists a neighborhood N s of [0, s) 
for which <f>(r, b) is analytic for r G N s . By the uniqueness of solutions to (|4ip these sets are 
nested, and in fact 0(r, b) is analytic in Af = U sg [ 0) _R b )A^ s . Assume Rb < 00, but 

limsup \<f>(r, b)\ < 00. 

Then, from the integral equation (j4"4"|) . it is easy to see that 0, 0' both have limits as r — )■ i?^. 
Using these limits as initial conditions for the differential equation, and applying Theorem 8.1 
of [CL] , there exists a C-neighborhood of Rb and an analytic solution of the equation, which 
extends <p(r, b). This contradicts the definition of Rb as the supremum. By restricting b and 
r to R, the complex analytic solution is real analytic on the desired domains. □ 
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We next prove some additional properties of F(r, R) 
restrict to r G 1. 



(r, bpt). In the following we 



Lemma 5.12. Suppose F(r) is a solution of ( l4"Tj) with < F(r) < 1 on (0,r ). Then: 



F"(0) = 0; 



F(r) 



< 0; 



(46) 



Proof. Since F(r)/r — > b as r — > 0, by the equation and L'Hopital's rule, 

rF'(r) - F(r) 1 



lim F"(r) = lim 

r->0+ r->0+ 



lim F"(r) 



Hence, F"(0) = 0. To verify the second conclusion, let h{r) = F(r)/r, and calculate 
1 (r 3 ti(r))' = h"(r) + -h'ir) = -F(r)(F 2 (r) - 1) < 0. 

In particular, r 3 h'(r) is strictly decreasing. Since r 3 h'(r) | r= o = 0, we have h'(r) = (F(r)/r)' < 
0, as claimed. 

Finally, for the third statement, 



F' 2F 



F\' fh\' h! 



h 



< 0, 



by the second conclusion. 



□ 



As the boundary-value problem ( )40|) does admit a unique solution for each R > (ob- 
tained as an energy minimizer,) for each R > 0, there is a unique choice of b = > for 
which <p(R, bji) = 1. The following facts are well-known, but we sketch the proof here: 

Lemma 5.13. Let 0(r, b R ) solve fl40]). Then b R > 0, < 0(r, 6 R ) < 1 and d r (j)(r, b R ) > for 
r E (0,R). 

Proof. The existence of a solution of ( HUj) follows from minimizing the energy, 



r r 



-(F'(r)) 2 + ^ 2 (r) + -(F 2 (r)-l) i 



r <ir, 



over F G H l {(0, R);r dr) with F(0) = 0, = 1. Since E{\F\) = E(F), \F\ is a 

minimizer if F is. By the strong maximum principle, we find that |i^|(r) > 0, and thus 
F > in (0, R\. Applying the maximum principle to the equation satisfied by F 2 (r) — 1 we 
may also conclude that F(r) < 1 for r 6 {0,R). We next prove that this F is monotone. 
First, we claim b = F'(0) > 0. Indeed, if F'(0) = 0, then F is a fixed point of the equation 
with 6 = 0. Since the solution is unique for each fixed b, we must have F(r) = 0, 
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which is a contradiction. Thus, b > 0. Suppose that F is not monotone. Then, there exist 
< R± < i?2 < R such that F has a local maximum at R\ and a local minimum at R2, with 
F(.Ri) > F(R 2 ). Therefore, we have F'{Rx) = = F'(i2 2 ), and F"(Ri) < < F"(i2 2 ), and 
so 

F ™ < F(J2 X ) (1 - F 2 (^)) , and ^ > F(i? 2 ) (l - F 2 (i? 2 )) . 



Ri R 



2 

Consequently, we have 

Rf < 1 - F 2 (i?i) < 1 - F 2 (i? 2 ) < R- 2 , 

a contradiction. Hence, F(r) is monotone, and F'(r) = d r <f)(r,bn) > for all r 6 (0,i?). 
Finally, uniqueness of the solution of fj40|) is proved in [M2j . 

□ 

We now consider the linearized operator, 

Lg(r) := -g" - X -g' + ^g - (1 - 3F 2 (r, i?))^, 
around F(r, R) = 0(r, 6^), the solution of (140]) . 

Lemma 5.14. Let Ai denote the smallest eigenvalue of L with Dirichlet boundary conditions 
on (0,R). Then X 1 > 0. 

Proof. Let g G ((0, -R)), g{r) > 0, be an eigenfunction associated to the first eigenvalue 
Ai, Lg = Xig. Now, h := F'{r) > solves a similar equation, and in fact 

IF F' 
Lh= — - — >0, 

by Lemma 15.121 Multiplying the equation for g by rh, the equation for h by rg, and 
integrating by parts, we have 

f R f R f2F F'\ 
Ai / g(r) h{r) r dr = j (— j g{r) r dr > 0. 







In particular, Ai > 0. □ 
We may now examine the dependence of the solution <p(r, b) on the shooting parameter 

b: 



Lemma 5.15. Let b = bu, with (p(R, 6r) = 1. Then ^0(r, bn) > 0. 
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Proof. Let h(r) = db<f)(r,b). Then we have 

- h" - -ti + — h = (1 - 30 2 (r, b R ))h, h(0) = 0, h'(0) = 1. (47) 

Suppose that h vanishes somewhere in (0, R\. Let Rq G (0, R) be the smallest value of r > 
for which /i(-Ro) — 0. Since h'(0) = 1, we must have R > and /i(r) > in (0,-R ). 
Choosing g{r) > the eigenfunction associated to the smallest eigenvalue Ai of L, 

-9" - -9 + \g = (1 - 30 2 (r, b R ))g + A^, 

we multiply by rh and integrate by parts on (0, Rq) to obtain: 

rRo 

Ai / g hr dr = 0. 



This is impossible, as g, /i > in (0, Ro), and hence h{r) = db<f>(r, 6r) > in (0, R]. □ 
Lemma 5.16. or is analytic in R> 0, and d^bR < 0. 

Proof. Since b R is defined as the solution to the equation (p(R, b) = 1, and from Lemma 15. 151 
we have dj,4>(R, 6r) > 0, applying the analytic version of the Implicit Function Theorem (see 

[Nil, we conclude that both dnbf> = — ^tttt r~r < and the dependence of b R on R is 
real-analytic. □ 

Proposition 15.81 now follows trivially from Lemma 15.111 and Lemma 15.161 as f\(r) = 
F(y/\r, Va), and F(r, R) = 0(r, is the composition of analytic maps in a neighborhood 
of the positive real axis A > 0. 

We may now apply analytic perturbation theory to simple eigenvalues of the operator 
Ai\ associated to the quadratic form Q\. 

Lemma 5.17. Assume that, for some Xq > and n > 1, /^(Ao) is a simple eigenvalue of 
M.\ . Then: 

(a) There exist 5, n > such that for X in a (complex) neighborhood D,s(Ao), the operator 
M.\ has exactly one isolated simple eigenvalue fin (A) G D^(/in ''(Ao)). 

(b) There exists a normalized eigenvector (ao(-; A), a2(-; A)) of Ai\ with eigenvalue fin\X), 
each depending analytically on X & 0,5 (A ). 

We remark that, by Lemma [5T7] the ground state eigenvalue /i[ (A) is a simple eigenvalue 
of A4\ for each fixed A > 0, and so the conclusions of Lemma 15.171 apply to /^(A) in 
particular. 
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Proof. First, by Proposition 15.81 f\(r) is real-analytic in both r and A > 0, and thus may be 
extended to complex A as an analytic function of both r, A, for A G ©5 (Ao) for some 5q > 0. 

Next, we observe that, for A G D,5 (Ao), M.\ is an analytic family in the sense of Kato, 
which has compact resolvent for all A G ID^Ao) and is self-adjoint for A > 0. The conclusions 
(a) and (b) of the lemma then follow from the Kato-Rellich Theorem (Theorem XII. 8 of 



Monotonicity of simple eigenvalues 

In order to study the dependence on A of the eigenvalues of the linearized operator M.\(ao, a-z) 
(defined in ( 158]) ). we proceed as in the proof of Lemma WM we replace the dependence on A 
by a dependence on the domain (0, R), via a change of variables. In this way, we define the 
quadratic form 

Qi?(a ,a 2 ) = Q^\a ,a 2 ), where a (r) = a {rR), a 2 (r) = a 2 {rR), R = \/A. 

The associated operator is then denoted by Mr. We observe that the eigenvalues c{M.r) = 
{/t n (i?)} ne N (ordered by the min-max principle, and repeated by multiplicity) are related to 
the eigenvalues fin (A) via fin \\) = R 2 fi n {R). 

Proposition 5.18. Suppose £t n (Ro) is a simple eigenvalue of for some Rq > 0. Then 



Proof. We use the family 0(r, b) of solutions to the Cauchy problem (14T!) above, and recall 
that 0(r, b) is real analytic in both (r, b). 

Let Ro > be fixed, and r\ > given. We recall that for any R > there exists a unique 
b = 6r > such that the solution to the boundary- value problem (HUj) is F(r, R) = 0(r, 6 R ). 
Furthermore, we assert that: 



EH]). 



□ 



p n (Ro) < 0. 



there exists e > so that |6r — &r | < r\ whenever \R — Rq\ < e; 
there exists C > so that b'(R) < —C < whenever \R — Rq\ < e. 



(48) 
(49) 



Indeed, both follow from the conclusions of Lemma 15.161 and Proposition 15.81 
Next, we claim that there exists a constant C\ > so that 



</>Mk) >c ir , 



for all r G [0, R], R G (Rq -e,R + e). 



(50) 



From statement (1481) above, and the analyticity of (ft, there exists 5 > such that 



<t>'{r, b R ) > ^ for all r G (0, 5) and \R - R \ < e. 
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Since <f>(r, b R ) > for r > and for all R > 0, 

k :— min 6(r, b R ) > 0. 

|K-iJ |<6 

Let C\ := min j-^ 2 -, ^ j. Then, putting the previous two estimates together we conclude 
that F(r,R) = (f>(r, b R ) > Ci r, for all r G [0,-R] and |i? — Rq\ < e, and the claim (!50|) is 
established. 

The next step involves the derivative db<p(r, b r ) =: /i(r, R). By Lemma |5.15[ /i is analytic 
in both (r,R), and h(r,R) > for all r G (0,i?], 6 > 0. We recall that h solves OH]), and 
ti(0,R) = 1 for all R. 

Following exactly the same arguments used in proving ( 15 Op . we obtain that the existence 
of a constant C<i > so that 

h(r,R)=db<f>(r,b R )>C2r, for all r e [0, R), R e (Rq - e, Rq + e). (51) 

As a consequence of (|5T|) and fT49|) . we have 

dnF(r, R) = d R <j)(r, b R ) = <%0(r, 6 fl ) 3 fl 6 B < -C C 2 r, 

for all r G [0, R] and |i2 — Rq\ < e, using statement (|49j) above. By the mean-value theorem, 
for any R G (Rq, Rq, + e) and for all r G [0, Rq], there exists R G (i?o, R) with 

F 2 (r, .Rq) - F 2 (r, i2) = -2(i2 - R )F(r, R)d R F(r } R) 

> (i?-i? )C 3 r 2 , (52) 

with constant C 3 = 2C CiC 2 > 0. 

We are now ready to bound the eigenvalue from below. Let R G (Rq, Ro + e) be fixed. 
By Lemma 15 .171 there exists 5 > for which p, n (R) is simple for \R — R \ < 6, and jx n (R) 
is analytic in that interval. Let E® denote the linear span of the first n eigenfunctions, 
Wi = (aj, a\), . . . , w n = (oq,^), of A4 Ro , each normalized with ||wfc||z, 2 = 1, k = 1, . . . , n. 
As each Wk(Ro) = 0, extending their definition by zero for r > R , each lies in the domain 
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1 

max 



of the the operator A4r for all R > Rq. Thus, by the Courant-Fischer min-max principle, 
p, n {R) = inf max Qr(w) 

dim E=n w eE 

IMl2 =1 

< max Qr(w) 

rtt 

Q Ro (w)+ / [F 2 (r,R)-F 2 (r,R )] [(a 2 + aj) + /3(a - a 2 f] r dr 

I™II2 =1 

< jl n {R ) - C 3 {R - Ro) min / (ao^i^ + a^r,^) 2 )?- 3 ^ 

IIHl2 = 1 

< frn(Ro) — C 4 (R — i? ), 

with constant C4 > independent of i£, using ( 152]) and the finite dimensionality of £7°. Since 
/i(-R) is isolated and simple in a neighborhood of Ro, by Kato-Rellich it is differentiable at 
Ro. By the above estimate, we conclude that fi' n (R ) < — C 4 < 0. 

□ 

From Lemma [5.71 we thus have: 

Corollary 5.19. Denote by fti(R) the smallest eigenvalue of Qr. Then ^[(R) < for all 
R>0. 

In order to return to the problem on a fixed ball D 1; with parameter A = R 2 , we recall 
that the eigenvalues of M.\ and M. R are related via ^n\\) = R 2 jl n (R), A = R 2 . So, 
j^fin\\) = ^Rfi' n (R) + fi n (R) is negative at a simple eigenvalue at the point at which 
fi n (R) = (that is, exactly at a bifurcation point): 

Corollary 5.20. Suppose fin\X) is a simple eigenvalue of M.\ for |A — A | < d~, and 
l£\\o) = 0. Then, ^fi { n\X ) < 0. 



Finally, we show that at least one eigenvalue of Ai\ must cross through zero as A increases: 
the ground state (A). 

Lemma 5.21. There exists a unique \p > so that /j,i\\) > for A < \ p and ^\\) < 
for A > X p . 

Proof. By Corollary 15 . 1 91 above . fi\{R) is strictly decreasing. From Corollary 15.201 it suffices 
to show that fti(R) < for some sufficiently large R. Thus fii^Rp) = at a unique Rp > 0, 



whence ^(Xp) crosses through zero at a unique A^ = R\. To do this we argue as i 



111 
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Theorem 2 of [Ml]. As R -> oo, the radial profile F(-,R) F^-) in C k ([0,R]) for all 
R > and /c G N, with i 7 ^ the modulus of the unique entire equivariant solution of the form 
Uoo = F OC) (r)e ie '. We have already shown that for < < 1, the entire equivariant solution 
= ^(ttoojttoo) is not a local minimizer (in the sense of de Giorgi) in M. 2 , so there exists 
R > and $ G Cg°(B B ) for which E^U^ + $;D fl ) < By an approximation 

argument, we could then conclude that £^({7r)[<I>] < for some $ G C£°(0#), and hence 
ti ( i\R 2 ) < for that value of R. 

Instead, we give a more direct proof, using Q R . Let d (r),d 2 (r) be the ground-state 
eigenfunctions. Define A = do — d 2 and S = do + d 2 . In terms of A, B we have: 

Q fl (d , d 2 ) = 4vr jf* j(A') 2 + {B'f + - A) 2 + (F 2 - 1)(A 2 + B 2 ) + 2/M 2 } r dr 
=■■ Qr(A,B). 

We would like to make the choice A = F^ and B = F^/r, but these functions are not 
admissible as test functions, since they do not vanish at r — R. Nevertheless, since F^ 
vanishes linearly at r = 0, A(r) and B{r) are regular near r = 0, and since A — B = r(F 0O /r)', 
the second term in Qr is well-defined. Moreover, it is well-known (see |HHj ) that the 
derivatives of F^ decay sufficiently rapidly as r — > 00 in order to have A and B in the 
domain of definition of Qoo, and 

Qoo(AS) = lim Q R (A,B). 

R—)-oq 

This last quantity we can evaluate exactly, using the equations which A, B solve (see (18') 
of [Ml] ): 

_ A " _ \ A < + %{A - B) - (1 - 3i£)A = 
- B" - -B' - - 2 {A - B) - (1 - Fl)B = 

ry y 

We multiply the first equation by A, the second by B, add, and integrate by parts to obtain 

POO 

Q oc (A,B) = 8n(P-l) / F^A 2 rdr<0, 

Jo 

when < /3 < 1. By approximation, we may find a large i? and i,Be C°°(0, i?) such that 
Q fl (i, B) = Qoo(A, B) < 0. Thus, for sufficiently large R, ti[ 1] (R 2 ) = R 2 fii(R) < 0. □ 

Bifurcation at simple eigenvalues 

We are finally ready to prove bifurcation of symmetric solutions to fl34|) at simple eigenvalues 
of L\. The following lemma summarizes the previous results on the eigenvalues of L\, 
obtained in the previous parts: 
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Proposition 5.22. Assume < (3 < 1. 

(a) Suppose /i„(A ) is a simple eigenvalue of M\ Q - Then: 

(i) /U n (A) is an eigenvalue of L\, with eigenspace 

X x = {sWf. seR, ^gS 1 }, 

with = (w£, —w^), vj£ = Wi = £ ao(r; A, n) — £ a2(r; A, n) e 2l6> . 

fii) 5oi/i the eigenvalue /i„(A) and normalized eigenvectors are analytic in a (complex) 
neighborhood of X . 

(iii) J//i n (A ) = 0, then /<(A ) < 0. 

(b) There exists a unique \p > /or which the ground state eigenvalue is given by /xi(A J g) = 
= Hi (\p), and is a simple eigenvalue of M.\„. 

Proof. Statement (%) follows from the reductions in Lemma [5.31 and ( 157)) . together with the 
description of the eigenspaces of C x in Lemma 15.61 The analyticity claimed in (ii) was 
proven in Lemma 15.17) and (iii) follows from Corollary 15.201 Part (b ) puts together the 
results of Lemmas 15.71 and 15.211 with that of Corollary 15.51 □ 



We are ready to analyze the bifurcation of solutions at isolated simple eigenvalues of A4\ 
which cross zero at A > 0. We write f eHasf = U\+V with V = (v + , vJ) G H^(Di, C 2 ). 
Then, \1/ solves (154")) if and only if 



with 



= F(V, A) := L\V + \H(V, A) 



H + (V,\) = (\v + \ 2 + \v^)(v + + ^=) + 2 



L+V + XH+(V, A) 
L^V + \H-[v,\) 



ux 
V2 



v + + n_ ) v + 



+ /3(K| 2 -K| 2 ) (v + + ^=) + 2(3( 



Ux 



V2' 



Ux 



H~(V, A) = (|n+| 2 + |«_r)(«_ + -£=) + 2 



ma 



\V2 

t> + + n_ ) n 



, f + — n_ ) n 



n_ > v 



(53) 



(54) 



v/2 \V2' 
-/3(K| 2 -K| 2 ) (n_ + ^)-2/3^,n + 

The above defines a smooth map F : #<J(Oi;C 2 ) x M+ -»■ ^(D^C 2 ). If we iden- 
tify a complex vector V = (n + ,n_) = (ni + in 2 , ni + in 2 ) G C 2 with the real vector 



,1 „,2 „,1 „,2 



(v+jf+j t>l,nij G M 4 , then by Proposition 15.81 we recognize that the symmetric solution Ux 
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is real-analytic as an M 4 - valued function of (x,X), and the map F is likewise real-analytic, 
viewed as a map of V G //o(Oi;IR 4 ) an d A. We also note that Hy(0, A) = 0, and so 
Fy{0, A) = L\, the linearization around the symmetric vortex solution U\. 

We now state our bifurcation result. Denote by X^ the orthogonal complement of X\ in 
^(D i; C 2 ). 

Theorem 5.23. Suppose Xq > is such that, for some n > 1, fin\Xo) is a simple eigenvalue 
of Ai\ with fin\Xo) = 0. Then (U\ , Xq) is a point of bifurcation for the equations (1341) . In 
particular: 

(1) there exists a neighborhood M of {U\ , Ao) in H x (0, oo), 5 > 0, and real analytic maps 
$ : (-8, 5) -> X^ o and : (-5, 5) -)■ (0, oo) with $(0) = ; 0(0) = A , such that there 
exists a non-equivariant solution (I34j) of the form 

= ^(t) + ^(^i + $(t)), A = <p(t) 

for all \t\ <5 and £ G S 1 . Moreover, #(-*,£) = i2_i*(t,0 and 0(-t) = 0(t). 

ylny solution A) o/ (1341) m the neighborhood M is either an equivariant solution 
(U\, A) or of the form (*&(£,£), 0(t)) above. 

(3) All solutions of (1341) wi A/" satisfy T^ = \^ (where the involution T is defined in (1351) .) 
Moreover, each component ip±(t,£) o/\l/(£, £) = (ip + (t,!;),ip-(t,(,)) has exactly one zero, 
and their zeros are antipodal and distinct from the origin. 

By (b) of Proposition 15.221 we may apply Theorem 15.231 at A = A/3, and obtain bifur- 
cation at the ground state eigenvalue /ii(A / g) = 0, which implies Theorem 11.51 stated in the 
Introduction. 

Proof. Given any V G Hq(JB>i;C 2 ), there exists a unique a > 0, £ G S 1 , and Z G 
with V = aW/: + Z = (aW 1 + Z), with Z = RjZ G and Wi = W ± (X) as in (%) of 
Proposition 15.221 normalized with ||Wi||£2 = 1. Define the Hilbert space y = {aWi} eg) X^ o , 
as a subspace of Hq(JB>i, C 2 ). We then consider the equation F(V, A) = as in ( J53l) . restricted 
to V G y. Since F is equivariant under the action, every solution V G y gives rise to an 
orbit of solutions in Hq(B>i] C 2 ), while any solution V G Hq(3i; C 2 ) corresponds to a solution 
R^V G J by an appropriate choice of £ G S 1 . Thus, it suffices to consider the equation (1531) 
in the smaller space 3^ to determine the solution space in Hq(B>i] C 2 ). 

We next list some properties of F restricted to y. First, it remains true that F(0, A) = 
for all A > 0, and F is a real analytic function of V G y, X > (by Proposition 15.221 
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(ii) and Proposition 15.81 thinking of U\,V G C 2 as real vectors in M 4 .) By the restriction 
to y, we have Fy(Q, Aq) = L\ with kery(X^ ) = {aWi}, and thus dimkery(L.\ ) = 1 = 
codimRan(LA ). Lastly, we calculate the derivative /x^(Ao) in terms of the function F: as 
/i n (A) = (Wi(A), F y (0, A) Wi(A)), by (raj of Proposition EM we have 

> // n (A ) = 2(W 1 '(A ), F v (0, Ao)W 1 (A )} + (Wi, Fy jA (0, A ) Wi) 

= (^,F V , A (0,A )W4), (55) 

as F v (0,X )W 1 = L^Wi = (since // n (A ) = 0). We now claim that F V)A (0, A )Wi(A ) £ 
Rany (L Ao ). Indeed, assume the contrary, so there exists X G y with Fy(0,\ )X = 
Fv,x(0, Ao)H / i(Ao), and take the scalar product with Wi(Xq). We have 

> // n (Ao) = (Wi(Ao), F v (0, A )X) = 0, 

a contradiction. 

The celebrated Crandall-Rabinowitz bifurcation theorem (Theorem 1.7 of |CR] ) may 
then be applied to F in the space y. We note that since F is an analytic map, by invoking 
the analytic version of the Implicit Function theorem in the proof of |CR] the maps obtained 
will be real analytic. We conclude that there exists a neighborhood Af of (0, Ao) in y x (0, 00), 
5 > 0, and real analytic maps <fi as in the statement of the theorem, such that 

F-^O} n Af = jV, A) = (t[Wi(A„) + $(*)], 0(*)) : |t| < 5} U { (0, A) G A/j . 

Define the neig hborhood TV G i?£(Bi; C 2 ) as the union of the images of Af under the action 
of £ G S 1 . The characterization of the solution set in Af in statements (1) and (2) then 
follows. 

It remains to verify the symmetry results in (1) and (3). First, we note that if (V 1 , <f){t)) 
is a solution in Af of the above form, then 

V := R-iV 1 = (-t)(Wt - $(*)) e Af ay 

is also a non-equivariant solution with A = (f>(t). Therefore, there exists s, \s\ < 5 for which 
(V,(f>(t)) = (V s ,<j)(s)). Since W l JL A^, we have 

(s + t)W x = s$(s) - t$(t) = 0, 

and hence s = — t and $(s) = $(— £) = —$(£). We conclude that 0(— t) = 4>{t) and 
V - * = i?_iV rt , which finishes the proof of (%). 
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Next, we show TV* = V* for the involution T. As above, TV* is also a solution with 
the same A = <p(t), in the neighborhood J\f. We note that TW\ = Wi, and by following the 
same arguments as above, T$(t) = $(t). Since TU\ = U\, we conclude T\l/(t,£) = 
for the whole family of solutions. 

Finally, we consider the zero set of each component of £). Since \fr(t,£) — > U\ in 
C 2 as t — >■ 0, and the equivariant solutions U\ have exactly one non-degenerate zero (the 
origin) in each component, the same must be true for for \t\ sufficiently small. Since 

is fixed by the involution T, the zeros of the components ip-(t,0 must be 



antipodal. By Proposition ^. 2[ only the equivariant solution U\ vanishes in both components 
at the origin, so the zeros of ip± must be antipodal and distinct. This concludes the proof 
of Theorem 15.231 □ 

Remark 5.24. (a) Given that the solution curves are analytic, we may expand them around 
the bifurcation point Ao and (in principle) obtain further information about the direction 
and stability of the bifurcating solutions. For instance, we may calculate higher derivatives 
of A = 0(t), and obtain A' = <p'{0) = 0, and 



'(0) 



(L A ,V"(0),V"(0))-2A^ / \W X 



Since L\ p is positive definite on X^ o , the sign of the numerator is not clear a priori, so 
numerical approximation may be necessary to determine the details of the bifurcation at A^. 

(b) In a similar vein, if we compute the quantity in ( 155]) directly, we obtain 
(F v>x (0, X )W, W) = An I [d x (A(/ A 2 - 1)) [a 2 + aj] + d x (\f 2 x ) [a - a 2 } 2 } r dr 

JO A=A 

Expressed in this form, it is not apparent whether this quantity is non-zero. Only by recog- 
nizing the connection to the derivative //„(A) are we able to apply the Crandall-Rabinowitz 
theorem. 

(c) While we know (from Lemma [5.2 II) that the ground state eigenvalue /ii(A) must cross 
zero for G (0,1), it is unclear whether any of the higher eigenvalues can lead to other 
bifurcations of the symmetric solutions. 

(d) A very general result by Rabinowitz |Ra] shows that bifurcation always occurs at eigen- 
values of any finite multiplicity in a variational problem. Although the form of equation 
assumed in [Ra] is somewhat different than our F(V, A), the result nevertheless holds true in 
our setting, although the conclusions of the bifurcation theorem are weaker than the state- 
ment obtained by using the simplicity of the eigenspace as in [CRj . In particular, one may 
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conclude that non-equivariant solutions exist in any neighborhood of U\ when = /i n (Ao) 
is a degenerate eigenvalue, and by analyticity (see [L]) the continua of solutions form finitely 
many analytic curves, but there is no complete characterization of the solution set as in (1 ), 
(2) of Theorem EM 
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